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Electromagnetic traps may be divided into two categories: those that trap neutral 
particles or atoms, and those that trap charged particles or ions. The trapping mech- 
anism for the former usually exploits the force experienced by a magnetic dipole in a 
nonuniform magnetic field. Typical configurations consist of a magnetic quadrupole 
field supplemented with a mechanism to reduce trapping losses near the field zero. 
Two well-established examples are the Ioffe-Pritchard trap [I], |2|, which has end coils 
ensuring a nonzero-field minimum, and the time-averaged orbiting-potential (TOP) 
trap 0, which employs a rotating magnetic field. Discussions of various neutral- 
particle magnetostatic traps are given in Ref. 0. Charged particles and ions can be 
trapped with the Paul trap ||, which involves an oscillating electric potential and 
can simultaneously trap particles of both polarities. An alternative is the Penning 
trap P, R], which combines electrostatic and magnetic fields and exists in many forms 



10]. A general discussion of ion traps is given in Ref. [ 1 1 



The quantum behavior of many particles confined in a trap has been the subject 



of much recent experimental and theoretical work |fL2|, [13]. Most investigations have 
centered on systems involving large numbers of bosons. In the present paper, we 
address instead some quantum-physics issues for certain traps containing one to a 
few hundred fermions, with the system in the ground state. We seek to describe the 
physics of an additional valence fermion added to such a system, perhaps in a highly 
excited state, using a relatively simple analytical model. Note that the focus of the 
present work is primarily theoretical issues involving quantum physics. In particular, 
we disregard experimentally important issues such as thermal effects flj], jl3| . 

Part of our interest in situations of this type stems from the analogy with Rydberg 
atoms, which have played a central role in the development of tools for understanding 
multi-electron atoms |L5j. Certain related developments could emerge in the context 
of multiparticle traps of the type we consider. Rydberg states of alkali-metal atoms 
are widely used for experimental and theoretical investigations, in part because the 
associated electronic core is relatively simple, forming a closed shell. For related 
reasons, in the present work we primarily consider systems of trapped fermions with 
the valence fermion lying outside a closed fermionic shell. The analogy with alkali- 
metal atoms is best for traps with a confining potential that is purely isotropic in 
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three dimensions. For certain values of the applied fields, the Ioffe-Pritchard and 
TOP traps provide examples that closely approximate this situation. 

Both the Ioffe-Pritchard and TOP traps capture magnetic dipoles ft aligned op- 
positely to an applied magnetic field B by drawing them into a region of minimum 
|B|. The potential energy is U(x,y,z) = /x|B|. Dipoles aligned with the field are 
expelled from the trap. 

The Ioffe-Pritchard trap consists of two coils with aligned symmetry axes in, say, 
the e z direction and four conducting bars equidistant from and parallel to the z axis. 
The coils carry parallel currents, while the bars carry alternately oriented currents. 
In cylindrical coordinates, the magnetic field to third order in p and z is 

B a 

B 1P (p,(f),z) = B c e z + B[p [cos 2<pe p -sin20e ] + -^ (z 2 - p 2 /2)e z - zpe p . (1) 

In this expression, the coils are taken to have radius R c , to carry currents J c e^, 
and to be positioned at z = ±A C . The bars are assumed to be located at p = Si 
with = 7r/4, 37r/4, 5tt/A, 7tt/4, and to carry currents of magnitude Ii. The other 
quantities are a characteristic length, area, field, and field gradient: l c = {A 2 -\-R 2 ) 1 / 2 , 
a c = 3(AA 2 — R 2 )/2, B c = p I c R 2 /ll, B[ = 2p h/^Sf. Trapped magnetic dipoles have 
potential energy 



U IP (p,(p, z) = p,B c 



+ 2\Bl li) P 



(2) 



If the condition B' 2 /B 2 = 3a c /l* is satisfied, as is readily feasible in practical situa- 
tions, then the potential Eq. (||) is isotropic: 

-2\ 



U lv {r) = iiBAl + T —\ , (3) 



where r 2 = l\ja c . Note that the orientation of the dipoles is anisotropic, for example, 
lying roughly along the — e z direction near the z axis. 

The TOP trap involves two parallel coils as above, but carrying oppositely oriented 
currents, and two further pairs of coils providing an additional rotating bias field. 
Let the quadrupole coils have radius R q , be located at z = ±A q , and carry currents 
=F/ g e^, respectively. Let the pairs of bias coils of radius Rb be located at x = ±Ab 
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on the x-axis and at y — ±A b on the y-axis, and let the corresponding currents be 
Ix{t) — hcosuibt and I y {t) = I b smuj b t. The time dependence creates a magnetic field 
vector that to lowest order lies parallel to the z = plane and rotates with frequency 
oj b . To second order in x, y, and z, the magnetic field is 



B TOP (x, y, z,t) = \B b cos w b t + B q x H -j- (3x — r ) cos uj b t — 2xy sin u;&£ 



/ B b CL b 

+ ( i?^ sin tUfct + + 4 (3y 2 — r 2 ) sin cu^t — 2xy cos cj^t 



B b a b 



- \2B' q z + [x cos u b t + y sin u b t] z \ e z , (4) 

where the characteristic length, area, field, and field gradient are l b = (A^ + i? 2 ) 1 / 2 , 
a b = 3(4A 2 - J R 2 )/2, B b = fi I b R 2 b /l 3 b , B' q = 3fi I q R 2 q A q /2(A 2 q + R 2 q fl 2 . The potential 
energy is obtained by expanding fi |B TOP | for small p, z and time averaging, assuming 
that the frequency is high enough so the dipoles do not move appreciably over one 
cycle. The time averaging removes linear and cubic terms, producing the trapping 
potential 



U TOP (p,z) = \xB b 




Ob_ 

21* 



z 2 



(5) 



to third order. If the physical parameters satisfy B' 2 /B 2 = 3a b /7l b , a readily attain- 
able condition, then the potential becomes isotropic: 

U TOP (r) = fiB b M + i J , (6) 

where r\ = 14l b /5a b . Note that the dipole orientation is anisotropic, being opposite 
to the time-averaged field Eq. (|J) and pointing away from the origin along the z axis 
and radially towards the origin in the z = plane. 

The Schrodinger equation for a particle in either isotropic trap is 

(-^V 2 + fiB + f iB r ^)ij(r,9A) = Eij(r,9A) , (7) 

where m is the dipole mass, Bq = B c or B b , and r s = r c or r b . Separating in 
spherical polar coordinates, we write i(}(r,9,4>) = (r /r)W(r)Y(9,4>), where Y(9,4>) 
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are spherical harmonics, r = (h/muj ) l l 2 , and luq = ftfiBo/mr 2 ) 1 ^ 2 . This choice puts 
the radial equation for W(r) into a convenient form: 

/ h 2 d 2 h 2 L(L+l) „ , „ ,\ Trr/ , . i 

where the angular momentum quantum number L takes values L = 0, 1, 2, . . . and the 
energy eigenvalues are £jv = /ii?o + huo(N + 3/2), with principal quantum number 
N = L, L + 2, L + 3, . . .. The wave functions are given in terms of generalized Laguerre 
polynomials: 

W N , L (r) = C N , L (r/r ) L+1 exp(-r 2 /2r 2 )L^i / L 2 } 2 (r> 2 ) , (9) 

where the quantities Cn,l are normalization constants. 

In what follows, it is useful to consider situations where a cloud of fermions forms 
a closed shell of energy levels in an isotropic trap. The number of fermions forming a 
closed shell below a particular value of N is determined by the degeneracies of states 
in the single-particle system. Examining the full wave functions and recalling that 
the allowed values of N increase in steps of two units, the degeneracy for a given N 
is found to be (N + 1)(N + 2)/2. Note that the doubling in atomic systems due to 
the two spin orientations has no analogue here because only dipoles oriented against 
the magnetic field are trapped. The number of particles completely filling levels less 
than or equal to iV is (N + 1){N + 2){N + 3)/6. 

In the present work, we seek to construct relatively simple analytical one-particle 
models for the valence fermion in cases with more than one trapped particle. In 
addition to their intrinsic interest, such models could be used to make analytical 
predictions of physical properties or could provide a favorable starting point for per- 
turbative and other calculations. The methodology applied here to generate effective 
models for various trap systems is related to that adopted in the development of a rel- 
atively simple analytical model for the valence electron in Rydberg atoms [Bl|. This 



model has been used in a variety of contexts |17], [18|, including recently the prediction 
of certain experimentally observable features of long-term revivals in Rydberg atoms. 

Consider a system of fermions caught in a trap, with one particle excited relative 
to the others. This valence fermion can be regarded as moving in an effective potential 
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created by a combination of the trapping fields and interactions with the other trapped 
particles. We include among these interactions the quantum effects from the Pauli 
principle, which prevents the valence fermion from occupying filled levels, and also 
interparticle forces that act to modulate the trapping fields. In constructing analytical 
models for the valence fermion, we address first the issue of incorporating effects 
from the Pauli principle and subsequently examine an analytical extension that could 
describe other interparticle forces. 

To illustrate the idea, consider the particular tower of states \N, L = 0, M = 0) 
accessible to a trapped valence fermion with angular quantum number L = and 
azimuthal quantum number M = 0. If only one particle is trapped, then N > 
with N even. If there is also a core of four fermions filling the energy levels below 
N = 2, then a valence fermion with L = is restricted by the Pauli principle to 
levels with even N > 2. Neglecting for the moment other interactions among the 
trapped particles, the valence fermion in each of the two situations can access states 
with identical eigenenergies, except that in the five-particle case the N = level is 
inaccessible. The issue of constructing an analytical effective one-particle theory to 
describe the L = states of the valence fermion in the five-particle case can therefore 
be rephrased as the problem of finding an analytically solvable effective trapping 
potential with L = energy eigenstates identical to those of the one-particle case but 
with the lowest one-particle level missing. 

Isospectral problems of this type can be treated in several ways. One approach 
might be the inverse method JL9| . In the present work, we choose instead an alterna- 
tive method with a definite physical interpretation that is both elegant and relatively 
simple, based on supersymmetric quantum mechanics with the superalgebra sqm(2) 
fl6| , \L% . This supersymmetric technique takes as input a Schrodinger hamiltonian 
H + for which the ground-state eigenenergy vanishes, and it determines via super- 
symmetry a complementary hamiltonian H~ . By construction, the eigenstates of H~ 
are degenerate with those of the original potential except for the ground state, which 
is absent. The eigenstates of the two systems H + and H~ are related by an explicit 
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map. The two associated Schrddinger equations can be written as 

R± ± = (_cP_ + fd£\ 2 ^ cPU_\ ± _ E± ^ ± 
2m \ dr 2 \dr J dr 2 J 

where V ;± are the eigenstates of H ± with eigenvalues E . The combinations of deriva- 
tives of U(r) generate associated potentials V ± (r). Since by assumption H + and hence 
V + (r) are given, the function U(r) can be found by solving a differential equation. 
The form of V~ (r) and H ~ can then be deduced. The existence of the degeneracy and 
explicit map between excited states in the two eigenspectra is a direct consequence 
of the super symmetry. 

This procedure can be applied to Eq. (g). Subtracting the ground-state energy 
El = fiB +hu)o(L+3/2) makes the lowest eigenvalue of the radial hamiltonian vanish, 
as required by the method. This subtraction depends on the angular momentum L, 
which therefore must be fixed. The input potential becomes 



^ + M = |^^ + ^oV-^(2L + 3) , (11! 



with eigenenergies E^ L = huj (N — L) and eigenfunctions W^ L (r) = W^lO")- The 
solution for the function U (r) is 

U(r)= 1 -(-Y ~(L + l)\n(-) , (12) 



2 \r J \r Q 
from which we obtain the form of H~: 
( h 2 d 2 h 2 (L + l)(L + 2) . „ „ hw 0/ A 

(13) 

The energy eigenvalues are E^ g L = hu (N s — L), where N s = L + 2, L + 4, . . ., and 
the eigenfunctions are W^ g L (r) = W / /v s _i i i + i(r). 

The derivation shows that for the L = case discussed above the eigenfunc- 
tions (r o /r)W^ g L = o(r)Y L =o,M=o(0,(p) = (r /r)W Na - lt L=i(r)Y L =o,M=o^^ 4>) are effec- 
tive one-particle eigenfunctions for the valence particle in a trap containing a total of 
five fermions. We emphasize that these are not conventional oscillator eigenfunctions, 
as can be seen from the index structure. 
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It might be tempting instead to model the valence fermion directly using con- 
ventional oscillator wave functions, based on a shell-model approach where the Pauli 
principle is incorporated by hand. However, the supersymmetric eigenf unctions have 
several advantages. Unlike the conventional oscillator case, for which the lowest state 
must be excluded by hand, the supersymmetric states form a complete set. Moreover, 
the lowest-state radial eigenfunction W^ g=2 L=Q (r) in our effective model has degree 
zero and hence zero nodes, as expected for the ground state of the five-fermion system. 
A conventional oscillator wave function would have one node instead. 

The potentials Vq~ and Vq differ by an inverse-square repulsive term: 

V -(r)-V + (r) = ^±- 2 +hu . (14) 

In the present context, the additional repulsion in plays the role of the Pauli 
principle by preventing the valence fermion from occupying the filled lower levels. 
Note that a change in the angular-momentum barrier would produce a similar effect 
on the potential but would not connect states with the same value of L. 

The above arguments for the case L = can also be applied to other values 
of L. For example, in the case where the valence fermion has L = 1 the method 
produces a relation between two effective one-particle models, the first involving two 
trapped fermions with one in the ground state and the second involving eleven trapped 
fermions with ten filling the levels below N = 3. 

Iterations of the method produce further relations. When the valence fermion 
has L = 0, for instance, the effective one-particle model for the system with five 
trapped fermions described by Eq. (|T^) can in turn be related to another model 
for a system with 21 trapped fermions, 20 of which fill the levels below N — 4. 
Implementing this mathematically requires another shift of the energy zero so that 
the ground state of the five-particle system has zero energy. Repeating the procedure 
produces a series of interrelations between effective one-particle models for the L = 
towers of states of systems with n = 1,5, 21, 57, . . . trapped fermions. A similar 
iteration for L = 1 generates connections between systems with n = 2, 11, 36, 85, . . . 
trapped fermions. The two sequences of numbers are generated by the formula n = 
1 + (N + 1)(N + 2)(JV + 3)/6, with N = —1, 1, 3, . . . or iV = 0, 2, 4, . . .. 
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The preceding discussion has largely disregarded effects of interactions between 
the trapped particles. If particle interactions are entirely neglected, the degeneracy 
of each fixed- N level means that the same effective one-particle model applies to the 
valence fermion in traps containing n& = d + (N + 1)(N + 2)(N + 3)/6 particles, 
where d = 1,2, ... , (N + 1)(N + 2)/2 is the number of fermions lying outside a 
closed shell. Although there may be special situations where particle interactions are 
relatively small, the possibility of additional interactions among the d valence fermions 
suggests that the best effective models would typically have d = 1. Similarly, the best 
effective models should also be ones describing small numbers of trapped particles and 
a relatively highly excited valence fermion. These considerations favor, for example, 
the model for the five-particle case with one highly excited fermion. 

Issues involving particle interactions could be addressed using a variety of standard 
methods, such as mean-field theory, perturbation theory, and Monte-Carlo methods. 
Given our present focus, we pursue here instead the possibility of incorporating in- 
teractions via analytical modifications to the effective models. 

Interactions typically shift the energy eigenvalues of the valence fermion. We treat 
this as a shift A = A(N,L) in the principal quantum number. Since the effective 
model produces apparent integer shifts in L, we define for convenience an integer I = 
I(L) and introduce the effective principal quantum number iV* = N + I — A. It is also 
convenient to define N s = N + 21, a quantum number analogous to the spectroscopic 
principal quantum number in atomic systems. We therefore have N* = N s — I — A, 
with new energy eigenvalues E^* = fiB + Kojq (N* + 3/2). 

The problem is to find a modification of the radial potential in the above effective 
one-particle models such that the eigensolutions of the corresponding radial equation 
remain analytical but are associated with the modified eigenenergies En*. This minor 
miracle can be accomplished by adding the effective potential 

^£(£+1^+1) 
2m r z 

to the operator on the left-hand side of Eq. (||), where L* = L + 1 — A. The resulting 
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differential equation, 

/ h 2 d 2 



+ ^ L * (L * r2 + ^ + t^o + iwmgr 2 ) W(r) = E N *W{r) , (16) 



\ 2m <ir 2 2m 

has eigenvalues E^* with iV* = L*, L* + 2, L* + 4, . . .. In terms of the functional 
form given in Eq. (|9|), the eigensolutions are Wn*,l*{t). Note that the special case 
A = reproduces all the complementary hamiltonians if ± discussed above when 
[iBq + hu>o(L* + 3/2 — 2S) is subtracted from both sides of Eq. (|I~6|). Thus, the radial 
equation for H + is recovered by selecting S = and 1 = 0, while that for H~ is 
recovered by selecting S = 1 and 1=1. The lath iteration of H + is obtained with 
S = and I = Iq, while the corresponding if - has S = 1 and J = Jq + 1. 

The analysis given here has some similarities to the derivation of a relatively 
simple analytical model for the valence electron in Rydberg atoms discussed in Refs. 
]16| , |I7|, |i~8|l . For example, our treatment of the L = states of isotropic trap systems 
with 1, 5, 21, 57, . . . fermions resembles that of the s-orbital states of alkali-metal 
atoms in the Rydberg case. Note that, since in the present trap context the valence 
fermion is typically neutral, its interactions involve a dipole moment rather than a 
monopole charge. This suggests that the effect of interactions may be significantly 
smaller than in the Rydberg case and the effective models correspondingly better. 
Note also that the form of the effective potential (|1^) is reminiscent of the special 
analytical model investigated in Ref. . 

In the atomic case, orthogonality of the model eigenfunctions arises because the 
quantum defects in, for example, alkali-metal atoms are asymptotically independent of 
the principal quantum number. Similarly, in the present context the eigenfunctions 
WN*,L*{ r ) form an orthonormalizable basis provided A is independent of iV with 
A < L+/+3/2. The N- independence implies that for each fixed L only one parameter 
is needed, corresponding to a simultaneous shift of the principal quantum numbers 
for that tower of states. Experience gained in the atomic case suggests that rapid 
asymptotic iV-independence is likely to suffice for applicability of the models here. 
However, it is unclear a priori which trap systems have this feature, as it depends on 
details of the many-body dynamics. Experimental investigations establishing energy 
spectra for the systems considered here would be of interest. The eigenfunctions 
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Wjf*,L*(f) could then be used to predict other quantities such as transition rates. 

In contrast to the Ioffe-Pritchard or TOP traps, which can be isotropic in three 
dimensions, the Paul and Penning traps are generically isotropic in two dimensions.^ 
Effective one-particle models can nonetheless be obtained via a similar approach. For 
brevity, we restrict ourselves to outlining the treatment of the case where interactions 
of the valence fermion are neglected. Effective one-particle models analogous to Eq. 
(|TE|) that allow for level shifts can also be constructed for the Paul and Penning traps. 

The Paul trap consists of a time-dependent potential given in cylindrical coordi- 
nates by (f>(p, <f), z, t) = 0, z) cosfipt, where <fi is 

0(P,0^) = ||(^-p72) (17) 

with characteristic voltage V p and length d p . A quantum solution for the one-particle 
case exists |T[| and can be investigated experimentally [p2 |, but for simplicity we 



consider here an alternative approach for large Q p ^> Q p = (\/2 \qV p \jmd 2 % ) 1 l 2 in 
which q<ft is approximated by a time-independent effective potential 



V = ^Lv<j ) -V<j ) = ±muj 2 p (p 2 + 4z 2 ) , (18) 

where m and q are the mass and charge of the trapped fermion and uo p = Q 2 /4Q P . The 
associated quantum problem separates with a wave function of the form \I/(p, <f), z) = 
(p p /p) 1 / 2 X(p)T(0, z), where p p = (h/mup) 1 ^ 2 . Introducing the quantum numbers 
K = 0, 1, 2, . . . and M = 0, ±1, ±2, . . . gives 

T M>K {<f>, z) = A K exp(*M0) exp{-z 2 /p 2 p )H K {V2 z/ Pp ) , (19) 

where the Hk are Hermite polynomials and the Ak are normalization coefficients. 
The two-dimensional radial equation is 



1 In a particular rotating frame, a single particle in a Penning trap with a special ratio of applied 
fields can experience an isotropic potential in three dimensions. The incorporation of rotation 
effects, including those on any fermion core, could then allow a treatment similar to that for the 
Ioffe-Pritchard and TOP traps. 
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The associated energy eigenvalues are E N>K = hu p (N+2K+2), where N = \M\ , |M| + 
2, \M\ + 4, . . ., and the eigenfunctions are 

X N , m ( P ) = C Nm (p/ Pp f Il+1/2 exp(V/2pJ)4j 2 'V|/ 2 (pVp?) , (21) 

where the Cn,\m\ are normalization coefficients. 

The Penning trap involves an electrostatic field of the form (|17D, along with a 
uniform magnetic field B = B p e z . Defining the axial frequency u z = ({qVpl/md 2 ) 1 ^ 2 , 
the cyclotron frequency uj c = \qB p \ /m, and = (cu 2 — 2u; 2 ) 1//2 , the one-particle 
hamiltonian for q > is 

h 2 

H = -t^V 2 + \mVL 2 p 2 + \mJlz 2 + \huj c id^ . (22) 

The equation separates via \I/(p, 0, z) = (p /p) 1 ^ 2 ^(p)©(6', z), where p — (h/mtUc) 1 ^ 2 . 
The p equation is 

/ h 2 d 2 h 2 M 2 - \ 



(23) 

where M = 0, ±1, ±2, . . . and K — 0, 1, 2, — The energy eigenvalues are E^,k,m = 
h(QN + 2uo z K - u; c M + VL + w,)/2, where iV = |M| , |M| + 2, |M| + 4, . . .. The full 
eigensolutions involve generalized Laguerre and Hermite polynomials: 

*n,kmM,*) = ^,|M|(p/po) |A/| exp[-|(^) 2 -i(^) 2 + i M0 

xL^l lMl/2 (kp 2 /2p 2 ) H K (z/z ) , (24) 

where k = Q/uj c , Zq = (h/mujz) 1 ^ 2 , and Cn,k,\m\ are normalization coefficients. 

Equations ( |2Ti| ) and ( |2"B| ) are both radial equations for a two-dimensional oscillator, 
and the same approach to effective one-particle models applies to each. For brevity, 
we treat primarily the Paul case in what follows. Subtracting the ground-state energy, 
the equation analogous to (|TT|) is 

Vfa = ^ ' M|2 p ; + \™y - ^ P (|M| + 1) . (25) 

The partner potential is Vj^-i = Vn{}| +1 + 2hu p . The difference Vj^i — Vj^i again acts 
as an additional repulsion that can be regarded as preventing the valence fermion 
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from accessing a filled lower level. The corresponding partial eigenfunctions are 
X+ m (p) = X Nm (p) and X- stlM] (p) = X Ns - U m\+i(p), where N s = \M\+2, \M\ + 4, 
. . .. The eigenspectra E~^i M , = hu p (N — \M\) and E^ a , M , = hu p (N s — \M\) are de- 
generate, except for the ground state. Similar expressions arise for the Penning trap. 

For both types of trap, the potential Vr^i depends on M and the subtracted 
ground-state energy depends on K (see Eqs. ( PUD and (£2^)). The effective one-particle 
models therefore describe towers of states \N, M, K) with fixed K and M. For the 
Paul trap, the number of spin- 1/2 states with energy less than or equal to E^,k is 

p f E(E + 2){2E - 1)/12 , E even , 
n ^-\ (^-l)(2£ + 3)/12, Eodd, (26) 

where E = E N)K jhuj p . The Paul-trap systems related by effective one-particle models 
are therefore those with 1, 3, 7, 15, 27, 45, 69, . . . trapped fermions. 

Similar considerations apply for the Penning trap. Note, however, that in this case 
the number of spin-1/2 states with energy less than or equal to E N>M ^ K depends on 
the frequency tuning. The magnetron motion is also unstable, so the corresponding 
quantum-number combination (N + M)/2 cannot be too large. 

In closing, we remark that the radial equations for all the traps considered in 
this paper can be mapped into various radial equations for Coulomb-type potentials 
24J, among which is the usual radial Coulomb equation in three dimensions. This 



suggests that under suitable circumstances the analogies drawn above between atoms 
and trapped systems might be more than heuristic. It would be amusing, for example, 
to identify a trap system for which the oscillator energy shifts correspond via a direct 
map to known quantum defects for an alkali-metal atom. 

We thank R. Bluhm, G. Gabrielse, and R. Pollock for useful comments. This work 
was supported in part by the United States Department of Energy under grant no. 
DE-FG02-91ER40661. 
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